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I.  Introduction 

Perturbed,  Linear,  tine- invariant  systems  of  the  fora 
x(t)  »  A (£)  x (t)  +  B(£)u(t) 


(1.1) 


are  the  focus  of  this  paper.  Here  x  and  u  are  n-  and  m-dimensional  state 
and  control  vectors  respectively;  £  is  a  small  positive  perturbation 
parameter;  and  Ate),  Bte)  are  analytic  at  £  *  0. 

We  shall  show  that  an  algebraic  approach  to  the  study  of  the  above 
system,  with  A(£) ,  B(£)  considered  as  matrices  over  the  (local)  ring  W 
of  functions  of  £  that  are  analytic  at  £  *  0,  leads  to  new  perspectives 
and  results  on  multiple-time- scale  behavior  in  this  system.  Proofs  are 
omitted  here,  but  relevant  ones  will  appear  in  the  final  version. 

We  begin  with  the  undriven  situation,  where 

x  •  A(e) x  .  (1.2) 

We  assume,  with  no  essential  loss  of  generality,  that  Ate)  is  nonsingular 
for  £  in  (0,  eQ),  for  some  eq  >  0;  the  situation  of  interest  is  where  A (0) 
is  singular.  Section  2  reviews  a  familiar  special  case  of  this,  namely 
the  class  of  two-time-scale  systems  extensively  studied  by  Kokotovic  and 
co-workers,  [1],  [2],  in  the  so-called  explicit  form 


(The  forms  actually  used  in  Cl] #  (2]  cam  be  brought  to  (1.3)  by  a  simple 
change  of  time-scale.)  The  section  then  outlines  the  extension  of  these 
results,  to  multiple- time-scale  decompositions  of  systems  of  the  form  (1.2) , 


that  have  recently  been  obtained  by  Coderch  et_  al. ,  [3]  . 

The  procedure  suggested  in  (31  for  extracting  and  displaying  the 
multiple-time-scale  structure  of  (1-2)  is  described  in  terms  of 
operations  such  as  projection  and  pseudo-inversion.  We  complement  that 
viewpoint  in  this  paper  by  presenting  a  slightly  more  concrete  version  of 
it  that  turns  out  to  be  very  fruitful,  our  procedure  is  the  natural 
generalization  of  the  one  used  for  (1.3),  and-makes  clear  the  role  of 
the  invariant  factors  of  Ate)  (in  the  ring  W) ,  a  role  that  is  suggested 
(but  not  developed)  in  [3].  This  algorithm  is. presented  in  Section  3;  the 
basis  for  it  is  a  transformation  of  Ate)  to  Smith  form,  which  is  outlined 
in  that  section. 

Section  4  turns  to  questions  of  feedback  control  of  time-scales  in 
the  system  (1.1),  assuming  state  feedback  of  the  form 

u  *  K(e)x  ,  (1.4) 

with  K(e)  again  a  matrix  over  W.  Noting  the  above  interpretation  of 
invariant  factors  of  the  system  matrix,  the  question  of  invariant-factor 
assignment  in  Ate)  +  Bte)K(e)  is  raised  and,  for  the  case  of  left-coprime 
A(e)  and  B(e),  answered  rather  completely.  To  actually  assign  time-scales 
by  such  feedback  requires  that  certain  stability  conditions  be  also  satisfied, 
but  our  invariant  factor  results  show  what  the  limits  are. 
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2.  3ackqround 

A  special  case  of  (1.2)  that  has  received  a  great  deal  of  attention 
in  the  control  literature,  see  [1] ,  [2],  is  the  system  (1.3).  It  is 
known  that,  if  and  A^  "  A^  ^  Ai2  are  nons^ni?u^ar'  the 
eigenvalues  of  (1.3)  occur  in  two  groups,  one  being  of  order  1  and  lying 
"close"  to  the  eigenvalues  of  A^,  and  the  other  being  of  order  £  and  close 
to  the  eigenvalues  of  £  (A22  -  A^  A^  ^  A^)  •  If  both  the  latter  matrices 
are  Hurwitz ,  then  the  system  exhibits  well-behaved  two-time-scale 
structure,  in  the  following  sense: 


[V'>- 

— 

x.,(t)  +  x.  (£t)  +  0(e) 

lx  IS 

a 

r— 

rt 

i 

x2s(et)  +  o(e) 

,  t  >  0, 


(2.1) 


where 


*lf  "  \l  Xlf  '  Xlf  =  Xl(0) 


xis  =  ^n"1  Ai2  x2s  '  and 


X2s  *  (A22  ~  A21  A11  A12}  X2s '  X2s(0)  =  X2(0) 


The  subscripts  s  and  f  denote  slow  and  fast  subsystems. 

For  some  appropriate  constant,  nonsingular  matrix  T,  it  can  be 
shown  that 


lim  sup  ilTe^2^  T  1  -  e  d  (j  *  0,  (2.2) 

£40  t>0 


where  the  block  diagonal  matrix  A^  is  given  by 
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A 


d 


A22  ‘  A21  A11  A12 


(2.3) 


(2.2)  and  (2.3)  provide  an  alternative  definition  of  what  it  means  to 
have  well-behaved  two-time-scale  structure. 

The  above  decomposition  has  found  significant  applications.  Two 

criticisms  that  may,  however,  be  noted  are:  firstly,  that  the  system  is 

assumed  given  in  the  explicit  form  (1.3);  and  secondly,  that  the 

nonsingularity  assumptions  on  A  ^  and  A^  restrict  the  system  behavior 

to  two  time  scales.  Recent  work  of  Coderch  et  al_. ,  [3]  ,  has  attempted 

to  address  both  objections.  It  starts  with  the  more  general  form  (1.2), 

assuming  however  that  a  Taylor  series  for  A(£)  is  given.  It  is  then 

shown  that  an  expression  of  the  form  (2.2)  holds,  with  a  different  T  and 

with  A.  of  the  form 
a 


e1^ 


m-* 


(2.4) 


The  above  govern  behavior  at  the  various  time-scales .  They  are  obtained 
through  a  rather  elaborate,  though  systematic  sequence  of  operations  on  the 
Taylor  series  coefficients  of  A(£),  involving  cascaded  projections  onto 
progressively  "slower"  subspaces.  The  convergence  in  (2.2)  is  proved  under 
a  so-called  semi-stability  condition  on  matrices  derived  from  the  Taylor 
series;  this  condition  implies  that  the  are  Hurwitz. 


With  this  as  background,  the  role  of  Section  3  may  be  stated  more 
clearly.  We  show  in  that  section  that  the  Smith  decomposition  of  Ate) 
over  the  ring  W  makes  possible  a  change  of  variables  in  (1.2)  that 
brings  it  to  what  can  be  termed  the  explicit  form  of  (1.2) .  This  form 
is  the  natural  extension  of  that  in  (1.3)  to  the  case  of  multiple-time¬ 
scales.  It  is  then  shown  that  a  simple  nested  iteration  of  the  familiar 
procedure  of  [1]  used  for  (1.3)  will  result  in  (2.2),  and  directly  give 
the  jL  of  (2.3);  it  is  assumed,  for  (2.2),  that  the  A^  are  Hurwitz.  We 
also  show  that  the  same  result  is  obtained  by  application  of  the  procedure 
in  [31  to  the  explicit  form  of  (1.2). 

While  the  above  approach  provides  some  valuable  additional  perspectives 
on  known  results,  the  real  pay-off  appears  in  the  results  on  feedback 
control,  described  in  Section  4, 


3.  An  Algebraic  Transformation  for  Multiple-Time-Scale  Decomposition 

A  transformation  of  (1.2)  that  more  explicitly  displays  its 

amenability  to  multiple-time-scale  decomposition  is  obtained  by  employing 

the  Smith  decomposition  of  A(£)  over  the  ring  W  of  functions  of  £  that 

are  analytic  at  0;  see  [4]  and  [5]  for  example. 

It  is  easily  seen  that  W  is  a  Euclidean  ring,  with  the  degree  of 

a  scalar  being  defined  as  the  order  of  the  first  nonzero  term  in  its 

2  3 

Taylor  expansion  (e.g.  £  +  £  +...  has  degree  2).  A(e)  therefore 

has  the  Smith  decomposition 

A(£)  -  P(£)D(£)Q(£)  ,  (3.1) 

where  P,  D,  Q  are  all  nxn  matrices  over  W;  p,  Q  are  unimodular ,  i.e. 
i  P (0)  1  0  and  ( Q<0)  j  ji  0;  and* 


where  0  <  i^  < •  •  •<  i^.  (We  have  used  the  assumption  that  A(s)  is  non¬ 
singular  in  the  neighborhood  of  0  in  writing  (3.2);  in  the  more  general  case, 
some  of  the  diagonal  terms  would  be  0) .  Actual  computation  of  such  decomposi¬ 
tions  is  discussed  in  [4]  and  [5] .  (In  the  terminology  of  [5] ,  what  is  required 
is  to  transform  A(£)  to  the  matrix  D(e)  Q(£),  which  is  "row-reduced  at  0", 
through  row  operations  embodied  in  P  ^(e.1.) 


*  The  identity  matrices  I  may  have  different  dimensions. 


We  term  (3.7)  the  reduced  explicit  form  of  (1.2) ,  since  it  is  obtained 
by  simplifying  the  explicit  form  (3.5).  (We  have  assumed,  with  no  loss  of 
generality,  that  i^  »  0;  this  can  always  be  obtained  by  a  change  of  time 
scale  in  (1,2)'.)  The  partitioning  indicated  in  (3.7)  will  be  explained 
shortly . 

The  rest  of  this  section  is  devoted  to  establishing  the  following: 

1.  The  system  (3.7)  can,  under  a  natural  set  of  stability  conditions, 
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be  decomposed  to  exhibit  well-behaved  multiple-time-scale 
structure.  The  decomposition  procedure  is  a  natural  extension 
of  the  familiar  one  of  [1]  for  (1.3) .  It  is  also  shown  to  yield 
an  equivalent  decomposition  to  that  obtained  by  the  procedure 
in  [3]. 

2.  When  (3.7)  has  a  well-behaved  multiple-time-scale  structure, 

(1.2)  has  this  same  structure  as  well,  and 


lim  sup  j|x  -  P(0)  z [ 
£i0  t>0 


0  . 


(3.8) 


To  establish  these  results,  note  first  that  the  indicated  partitioning 
of  (3.7)  puts  it  in  the  form  (1.3),  so  that  the  familiar  decomposition 
procedure  of  [1],  [2]  can  be  used.  We  then  get,  as  in  (2.1) ,  the  following 
decomposition  into  slow  and  fast  subsystems: 


’z1(t) 

V,(t)  +  z,  (et)  +  0(e)" 

If  Is 

z .  (t) 

3 

z_  (et)  +  0(e) 

L  2,m 

L  2s  J 

(3.9a) 


where 


‘2,m 


LZmj 


(3.9b) 


Zlf  *  A11  Zlf  ' 


(3.9c) 


Z .  *  —  A,  .  [A,  ...  A,  ]  Z_  , 

Is  11  12  lm  2s 


( 3 . 9d) 
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’2s 


22 


m-1- 


m2 


2m 

V* 

£  A 


mmJ 


'2s 


and 


(3.9e) 


Ail 


( 3 . 9f) 


The  decomposition  holds  under  the  condition  that  the  subsystems  in  (3.9c) 
and  (3.9e)  are  stable. 

Observe  now  that  the  system  in  (3.9e)  is  itself  in  reduced  explicit 
and  may  be  subjected  to  the  same  procedure.  Iteration  of  this  (Schur 
complementation)  procedure  leads  to  a  decomposition  into  m  subsystems  of 
form 


*(i)  .  S  ,(i) 

Z  *  A^  Z 


i  »  1  to  a, 


(3.10) 


where  the  A^  are  given  by 
(i) 

*i  ■  Au  • 


with 


(JO 


13 


(k-l) 


-  A 


Oc-1) 

ik 


(k-l) 

*kk 


-1 


(k-l) 


and 


(0) 

A.  . 

13 


A. 


(1) 


The  A^  in  (3.9e)  are  actually  A^  in  the  present  notation,  and  z^f 
of  (3.9c)  is  z(1)  of  (3.10).  Also,  the  A^  in  (3.10)  are  precisely  those 
referred  to  in  (2.2),  (2.4).  Under  the  assumption  that  the  are  Hurwitz, 
we  find  that  (3.7)  has  well-behaved  time-scales.  The  number  of  variables 
of  (3.7)  at  each  time-scale  is  precisely  given  by  the  degrees  i^  of  the 


-m. 


invariant  factors  of  A(e)  (which  appear  on  the  diagonal  of  (3.2)). 


The  results  of  Section  3  have  established  the  role  of  the  invariant 
factors  (i.e.  entries  of  D(£)  in  (3.2))  of  the  matrix  A  (£ )  in  determining 
the  time-scales  of  the  undriven  system  (1.2).  "or  the  driven  system  (1.1), 
it  is  now  natural  to  ask  what  freedom  there  is  in  (re-)  assigning  time-scales 
by  application  of  the  state  feedback  of  (1.4) .  This  feedback  yields  the 
closed-loop  system 

X  *  F(£)x,  F(£)  *  A(£)  +  B  (£)*  K(£)  .  (4.1) 

A  key  question,  then,  is  the  question  of  invariant  factor  assignment: 
what  freedom  is  there  in  assigning  the  invariant  factors  of  F(£)  by  choice 
of  K(e) ?  The  following  theorem  provides  a  result  in  this  direction. 

Theorem:  Assume  that  A(£) ,  B(e)  are  left  coprime,  i.e.  that  [A(0)  3(0)] 
has  full  row  rank.  (Recall  that  A  and  B  have  dimensions  nxn  and  nxm  respectively.) 
Let  b  denote  the  rank  of  B(0) . 

1.  F(£)  can  have  no  more  than  b  non- unit  invariant  factors. 

jl  jb 

2.  There  exists  a  I<(£)  such  that  F(e)  has  £  ,...,£  as  its  non-unit 

invariant  factors,  for  arbitrary  non-negative  integers 

00 

(with  the  convention  that  £  =  0) . 

The  existence  of  well-behaved  multiple-time-scale  structure  in  (4.1)  that 
corresponds  to  the  above  invariant -factor  structure  can  be  guaranteed  if  the 
various  r\  (defined  as  the  were)  turn  out  to  be  Hurwitz. 

Some  results  are  also  available  for  the  case  of  non-coprime  A(e) ,  3(e) . 

In  this  case,  F(£)  is  of  the  form 

F  (£ )  *  W(£)F(£),  F(£)  *  A  (£ )  +  5(£)K(S)  ,  (4.2) 

where  W(e)  is  a  greatest  common  left  divisor  of  A(e) ,  3(£),  and  A(£) ,  5(e) 
are  left  coprime.  If  the  invariant  factors  of  F(e) ,  W(e)  and  F(£)  are  denoted 


-13- 


by  f .  (£)  ,  w^(e)  and  £Az) ,  and  ordered  such  that  the  i-th  one  divides 
the  (i  l)-th  one,  we  will  have 

wi(£)jfi(e)  and  5(e) ! f  (e)  -  (4.3) 

The  first  divisibility  condition  in  (4.3)  shows  that  every  invariant 
factor  of  F(E)  must  contain  the  corresponding  invariant  factor  of  W{£) . 

The  f^(E)  are  governed  by  the  above  Theorem,  and  conclusions  about  the 
f^(£)  can  then  be  drawn  from  the  second  divisibility  condition  in  (4.3). 
Conclusion 

A  promising  basis  for  an  algebraic  treatment  of  time-scale  structure 
and  assignment  in  linear,  time- invariant  systems  has  been  presented.  A 
wide  range  of  research  questions  has  thereby  been  exposed,  and  preliminary 
results  have  been  outlined. 
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